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Abstract 

Photoproduction at HERA is studied in ep collisions, with the ZEUS detector, for jp 
centre-of-mass energies ranging from 130-270 GeV. A sample of events with two high-py 
jets (pr > 6 GeV, r) < 1.6) and a third cluster in the approximate direction of the electron 
beam is isolated using a clustering algorithm. These events are mostly due to resolved 
photoproduction. The third cluster is identified as the photon remnant. Its properties, 
such as the transverse and longitudinal energy flows around the axis of the cluster, are 
consistent with those commonly attributed to jets, and in particular with those found for 
the two jets in these events. The mean value of the photon remnant pt with respect to 
the beam axis is measured to be 2.1 ± 0.2 GeV, which demonstrates substantial mean 
transverse momenta for the photon remnant. 
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1 Introduction 



Hard scattering between a real photon and a proton is expected, in lowest order QCD, to 
occur by two different mechanisms (Fig. 1). The photon may interact directly with a quark 
or gluon from the proton or it may resolve into its constituent quarks and gluons which then 
interact with the partons from the proton. These two processes are called direct and resolved 
photoproduction, respectively [[!]. 

In direct photoproduction, the entire photon participates in the hard scatter. In resolved pho- 
toproduction on the other hand, only a fraction of the momentum of the photon is involved in 
the hard scatter while the remaining momentum is carried by spectator partons. These par- 
tons fragment into a photon remnant which is expected to be approximately collinear with the 
original photon. The presence of resolved photon interactions at HERA has been demonstrated 
[0, |3| and the existence of the photon remnant has been confirmed. The separation between 
the direct and resolved contributions has also been reported, together with the measurement 
of a differential dijet cross section |3J]. Inclusive jet ^ [J and dijet J7J cross sections have given 
further information on the kinematic properties of these two processes and on the gluon content 
of the proton as well as on the parton structure of the photon. 

The presence of a photon remnant in photon-proton collisions has been used as a means to 
identify resolved photoproduction interactions. The study of the photon remnant itself, how- 
ever, is also interesting: a detailed comparison with leading-order (LO) predictions has not 
been performed and little is known about the internal structure of the photon remnant. On 
the theoretical side, although the point-like coupling of the photon to quark-antiquark pairs is 
included in parameterizations of the photon structure function, most Monte Carlo simulations 
model the resolved photon as a hadron, with collinear incoming partons. This results in a pho- 
ton remnant with low-pr with respect to the beam axis. Several studies have suggested that 
next-to-leading-order contributions or fluctuations of the photon into quark-antiquark pairs 
with high virtuality may lead to a 'photon remnant' which has sizable transverse momentum 
with respect to the incident photon direction || ||, [10 . 



In this paper, the photon remnant is isolated for the first time. Its properties are studied using 
a clustering algorithm and are found to be consistent with those commonly attributed to jets. 
These remnant jets are then compared with the jets emerging from the hard interaction and 
with LO Monte Carlo expectations. The data were collected with the ZEUS detector during 
the 1993 data-taking period. The study is conducted for -/p centre-of-mass energies (W 7P ) in 
the range 130 < W yp < 270 GeV. 



2 Experimental setup 



Details of the ZEUS detector have been described elsewhere The primary components 

used in this analysis are the calorimeter and the tracking detectors. The uranium-scintillator 
calorimeter [12| covers 99.7% of the total solid angle. It is subdivided into electromagnetic 



and hadronic sections with cell sizes of 5 x 20 cm 2 (10 x 20 cm 2 in the outgoing electron 
directionP]) and 20 x 20 cm 2 , respectively. It consists of three parts: the rear calorimeter 
(RCAL) covering the region —3.8 < i] < —0.75, the barrel calorimeter (BCAL) covering the 
region —0.75 < 77 < 1.1 and the forward calorimeter (FCAL) covering the region 1.1 < 77 < 4.3. 

x The Z axis is defined to lie along the proton direction; the Y axis points upward; the pseudorapidity, 
j] = — ln(tan |) where 8 is the angle with respect to the Z axis. 
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i(35) %/y/E{GeV) 
© 0.5 ns, where E 



The calorimeter has an energy resolution achieved in test beams of a/E = It 

for electrons (hadrons). The timing resolution for each cell is a t =1.5/ yE 
(GeV) is the energy deposited in the cell. 

The tracking system consists of a vertex detector [[H| and a central tracking chamber |14| inside 
a 1.43 T solenoidal magnetic field. The interaction vertex is measured with a resolution along 
(transverse to) the beam direction of 0.4 (0.1) cm. 

The luminosity is measured, using the electron-proton bremsstrahlung process, by electron and 
photon lead-scintillator calorimeters [15 installed inside the HERA tunnel. In 1993 the beam 
energies at HERA were E e = 26.7 GeV for the electrons and E p = 820 GeV for the protons. 
Typical electron and proton currents were about 10 mA and the instantaneous luminosity was 
about 6 x 10 29 cm~ 2 s . HERA operated with 84 colliding bunches. Additional electron and 
proton bunches circulated without colliding and are used for background measurements. 



3 Trigger and data selection 

In this analysis, photoproduction events are defined by requiring that the electron was scat- 
tered at small angles and was not detected in the calorimeter. This requirement corresponds 
approximately to a cut of Q 2 < 4 GeV 2 , giving a median Q 2 of about 10 -3 GeV 2 [Q. The 
trigger selects hard scattering events at low Q 2 . 

The ZEUS detector uses a three level trigger [II]. In the first level trigger, the calorimeter 



cell energies were combined to define regional and global sums which were required to exceed 



given thresholds [IB]. The second level trigger mainly rejected beam-gas interactions using 
timing information from the calorimeter. The third level trigger performed further rejection of 
beam-gas and cosmic ray events using information from both the calorimeter and the tracking 
chambers. An event was rejected if no vertex was found by the central tracking chambers or 
if the vertex was located in the region \Z\ > 75 cm. To reject beam-gas interactions, events 
were selected based on the following kinematic cuts 0: E tot — pz > 8 GeV, pz/E tot < 0.94 
and E^ ne > 12 GeV, where the calorimetric quantities E to t, pz and E^ ne are the total energy, 
the total longitudinal energy and the transverse energy excluding a cone of 10° in the forward 
direction, respectively. About 470,000 triggers were collected with these trigger conditions. 
As in previous studies of hard photoproduction M, |7J, the following offline cuts were applied to 
select the final event sample. 

• Beam-gas interactions were reduced by tightening calorimeter timing cuts, as well as cuts 
on the correlation between the vertex position (defined by the tracking chambers) and 
the calorimeter timing [ |i7|[ . 

• The E!j? ne cut was raised to E^ ne > 15 GeV to select hard scattering events. 

• To reduce beam-gas interactions, the event was rejected if less than 10% of the tracks 
pointed toward the vertex. 

• Deep inelastic scattering (DIS) neutral current events were removed from the sample as 
described in our previous publications [fj. f]J . 

• The fraction of the initial electron energy carried by the almost real photon, y = E y /E e 
where E 1 is the photon energy, was measured using the Jacquet-Blondel |1| estimator 
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of the Bjorken-T/: 

VjB 2E e 

The sum runs over calorimeter cells with energy E l and longitudinal energy p l z . To reduce 
uranium noise, the cell energies were required to be greater than 60 MeV (110 MeV) for 
the electromagnetic (hadronic) cells. This calculation assumes that the scattered electron 
was not detected in the calorimeter. For DIS events in which the electron deposits energy 
but is not identified in the calorimeter, y JB will be near unity. Therefore, for further 
rejection of DIS events, we required y JB < 0.7. To reject proton-gas interactions, y JB 
> 0.2 was required. These requirements correspond approximately to 0.2 < y < 0.85. 

• To remove charged current background and cosmic ray showers, a cut on ^t/V^t < 1-5 
GeV 1 / 2 was imposed, where j/y is the total transverse momentum and E? is the total 
transverse energy of the event. 

After these selection cuts, a sample of 99,894 events remained, corresponding to an integrated 
luminosity of 0.55 pb _1 . The estimated proton-gas and DIS background contributions are 0.4% 
and 1 to 2%, respectively. Cosmic ray and electron-gas backgrounds are negligible. 



4 Monte Carlo simulation 

In the following, the data are compared to Monte Carlo simulations based on the PYTHIA 
5.6 |19| event generator which includes leading-order QCD calculations. The HERWIG 5.7 [f20f 



event generator was used to check the PYTHIA results. The cut on the minimum transverse 
momentum, pTmin, of a hard scatter was set at 2.5 GeV. In PYTHIA, the photon flux is 
calculated using the Weizsacker- Williams approximation. The parton densities used were GRV 
LO |21j for the photon and MRSD_ [^] for the proton. For comparison, we also used the 
parameterization LAC1 |23j for the photon. 



The generated events were passed through a detector simulation based on GEANT 3.13 |24 
The same reconstruction program that was used in the data analysis was applied to the gener- 
ated events. The generated Monte Carlo event sample was obtained by combining the resolved 
and direct samples in proportion to the generated Monte Carlo cross sections (approximately 
7:1 for the GRV photon parton densities). 

In PYTHIA the distribution of the intrinsic transverse momentum, kt, of the partons in the 
proton and in the resolved photon is parameterized by the distribution, dN/dkf cx e~ kt ^ k °, 
where ko is a parameter which determines the hardness of the kt spectrum. The default value 
of ko for both the proton and the resolved photon is 0.44 GeV. An option in PYTHIA allows 
events to be generated using a different functional form for the k t spectrum and a different value 
of ko- This option has been used to generate events with a harder kt spectrum (see section 8). 



5 Analysis using the k± algorithm 

Previous analyses of photoproduction at HERA have implemented a cone algorithm to find 
jets 0, ^, |j. This algorithm [^] uses a cone of fixed radius in pseudorapidity and azimuthal 
angle space and maximizes the transverse energy within this cone. It is therefore well suited 
for high transverse energy jets. The photon remnant, however, is expected to deposit energy in 
the electron direction with low transverse energy. Due to the rapid variation of pseudorapidity 



3 



in this region, and to the low transverse energy, a cone algorithm cannot be used to identify 



the photon remnant. Therefore, we chose to use the k± clustering algorithm | 26fl . The analysis 
was done in the laboratory frame instead of the 'jp frame due to the uncertainties in the boost. 
The k± algorithm finds jets by iteratively merging clusters. Initially, clusters are individual 
calorimeter cells. For the Monte Carlo events, the algorithm is also used to cluster generated 
particles (see section 7). In the merging procedure, the quantity k± is evaluated for each pair 
of clusters, 

k ± = 2E^ in {L-coa6 nm ), 

where 8 nm is the angle between the two clusters n and m, and E m i n is the minimum value 
of the cluster energies E n and E m . In the small angle approximation, k± is the transverse 
momentum squared of the lower energy cluster with respect to the higher energy cluster. A 
pseudo-particle with infinite momentum along the Z axis is included in the clustering procedure 
to take the proton remnant into account. The value of k± between the pseudo-particle and the 
other clusters is calculated using the same formula as above. When all of the k± values have 
been calculated, the two clusters with the lowest kj_ value are merged. The four-momentum of 
the new cluster is the sum of the four-momenta of the two merged clusters. The calculation 
of k± is then repeated, replacing the two merged clusters with the new cluster. The iteration 
continues until the energy-angle resolution variable, Y (Y = kj_/Ej., where E? is the total event 
energy transverse to the beam axis), becomes larger than some threshold, Y cut . The value of 
Y cu t may be fixed or chosen on an event-by-event basis. For high transverse energy jets, the 
k± cluster-finding algorithm gives results which are similar to those obtained with the cone 
algorithm. 



6 Photon remnant identification 

To identify the photon remnant, we begin with an intuitive approach which is later justified 
using simulated events. As illustrated by Fig. lb, resolved hard photoproduction events have a 
final state which includes two high p? jets from the hard scatter as well as photon and proton 
remnants. Since we expect to find three clusters (in addition to the proton remnant) in each 
event, we choose the value of Y cut on an event-by-event basis so that three clusters are found in 
each event \Tf\\. These three clusters then should correspond to the two clusters from the hard 



scatter and the photon remnant. Since the photon remnant is expected to have low transverse 
momentum with respect to the beam axis, the separation between the photon remnant and the 
two jets from the hard scatter can be achieved, to a first approximation, by associating the 
photon remnant with the cluster having the smallest transverse momentum. For direct events, 
where we do not expect to see a photon remnant, the lowest transverse momentum cluster will 
either be part of the proton remnant or part of one of the two high transverse momentum jets. 
Therefore, the pseudorapidity distributions of the lowest transverse momentum clusters will be 
different for resolved and direct events. 

In Fig. 2a-c, we show the (uncorrected) pseudorapidity distributions of the three clusters ob- 
tained with the k± algorithm. The data (full circles) and Monte Carlo events (histogram) are 
shown normalized in the region r] cal < 1.6. While the two clusters with the highest p T cal (Fig. 
2a,b) are mostly found in the r] cal > region, the third cluster, with the lowest pt ccl1 , (Fig. 2c) 
is observed mostly in the rf al < region, i.e. in the photon direction. Data and Monte Carlo 
expectations agree for the two highest p T cal clusters except in the forward region, r] cal > 1.6, 
where we observe, in the data, an excess similar to that already reported in our previous anal- 
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yses [|, |7J. The peak observed in Fig. 2c, for the third cluster, in the negative pseudorapidity 
region is accounted for by the Monte Carlo simulation including resolved and direct processes 
(full histogram). The Monte Carlo distribution, however, is shifted slightly to lower values. 
The direct process, which does not contain a photon remnant, does not contribute to this peak 
as shown by the dashed line in the figure. Therefore, we are justified in using this method to 
separate direct and resolved events. In addition, the third cluster, when it is in the negative 
pseudorapidity region, can be associated with the photon remnant. 

In order to maximize the possibility that the two highest p T cal clusters stem from the partons 
in the hard scatter and to minimize the possibility that one of the two highest p T cal clusters 
is, in fact, part of the proton remnant, we require that the two highest p T cal clusters have 
high transverse momentum {prTl > 5 GeV) and that they are well separated from the forward 
region (rjf 1 ^ < 1.6). 

The distribution of the pseudorapidity of the third cluster, rj^ al , after the above cuts and the 
requirement E^ al > 2 GeV, is shown in Fig. 2d. The comparison with the distribution predicted 
for direct processes shows that the events with rj^ al < — 1 are almost exclusively due to resolved 
processes. The agreement of the Monte Carlo simulation for resolved plus direct contributions 
with the data is not perfect; in the data there are somewhat fewer events with large negative 
rjf- 1 values. The difference between the data and the Monte Carlo simulation is not improved 
when the photon parton parameterization LAC1 (dotted line) is used instead of GRV LO. For 
the following analysis, the resolved event sample is selected by requiring r]^ al < — 1. 
A measurement of the mean value of Y cut also demonstrates the difference between the two 
event samples on either side of rjf 11 = —1. Higher values of Y cut indicate energetic clusters 
which are spatially well separated. For the events at high r)™ 1 , the mean value of Y cut is 0.028. 
For the low rj^ al events, the mean value of Y cut is 0.063. This result suggests that the high rj^ al 
events contain a significant number of direct, two-jet events, for which the clustering procedure 
has been prematurely stopped. For the sample at rj^ al < — 1, however, the separation between 
the three clusters is quite distinct. 

A total of 1370 events satisfy the cuts. For these events, the fraction of the photon momentum 
involved in the hard scattering, x 7 , measured from the two highest p™ 1 jets H, peaks at low 
values as expected for resolved processes (not shown). By defining direct events, in this sample, 
as events with x 1 > 0.75 0, we obtain a contribution from direct processes of about 8%. The 
contribution from direct photon interactions, estimated using LO direct Monte Carlo events, is 
also 8%. 

7 Correction procedure 

In the following, the data are corrected back to the hadron level using the PYTHIA Monte 
Carlo program previously described. For the Monte Carlo events, the k± algorithm is applied 
independently at both the generated hadron level and the calorimeter cell level. In both cases 
the resulting clusters are sorted according to px- The detector level cuts in r] cal , pT cal , E^ al , 
and djb correspond approximately to two hadron jets with pri,2 > 6.0 GeV and 771,2 < 1.6, a 
remnant cluster with 773 < —1 and E$ > 2 GeV, and 0.2 < y < 0.85. Therefore, the hadron 
level Monte Carlo distributions were determined for events within this kinematic region. 
The correspondence between hadron and calorimeter clusters was determined by comparing 
the value of k± for each pair of hadron and calorimeter clusters. Each hadron cluster was then 
matched with that calorimeter cluster with which it had the lowest value of k±. After all cuts 
were applied at both the hadron and calorimeter levels, and neglecting the possible interchange 
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of the labels, cluster 1 and cluster 2, due to the sorting (this occurs in events in which the 
two highest px jets have similar p?), all three clusters were found to be correctly matched for 
97% of the events. For 3% of the events, a high-p^ cluster at one level was associated with the 
proton remnant at the other level. In less than 0.5% of the events, the lowest pt clusters at the 
hadron and calorimeter levels were not matched. 

The experimental shifts and resolutions of the measured variables were evaluated using the 
matched hadron and calorimeter clusters. The average measured value of rj is shifted, with 
respect to the hadron level, by less than — 0.02(+0.08) and has a resolution of 0.07(0.43) units 
for the first two jets (third cluster). The measured transverse jet momentum, with respect to 
the beam, is reduced, on average, by 15% and has a resolution of 11%. This is due either to the 
magnetic field acting on low energy particles or to dead material in the apparatus. The value 
of y is reconstructed with an average shift of —0.14 units and a resolution of 10%, also due to 
losses in the beam pipe and detector effects. 

On average, the third cluster contains 75 ± 20% of the photon remnant energy, as defined by 
the third cluster energy at the hadron level. Inactive material in front of the calorimeter results 
in an energy loss of about 20%, independent of the cluster energy. Particles lost in the beam 
pipe account for the rest of the energy loss; this effect increases with cluster energy, becoming 
comparable to the detector effects at measured energies above 10 GeV. 

The data were corrected with the following procedure. First, the contamination from events 
outside the kinematic range was estimated using Monte Carlo events and was subtracted bin- 
by-bin from the measured distributions. The resulting distributions were then corrected with a 
correlation matrix that was generated using the matched hadron and calorimeter clusters. The 
unfolding procedure is described in [2J] and includes acceptance corrections for the trigger and 
cuts described in section 3. The statistical errors were estimated by randomly varying, within 
their statistical errors, both the experimental data (before the background subtraction) and 
the Monte Carlo correlation matrix and calculating the root mean squared deviation in each 
bin. 

The systematic uncertainties were estimated by considering the following effects: different 
photon parton densities, a 5% energy scale uncertainty in the Monte Carlo simulation of the 
energy response of the calorimeter, a variation of the cuts on the measured quantities (y JB , 
PtTL Vi a L V3 11 )} corrections using the HERWIG Monte Carlo simulation, an increase in the 
direct cross section by a factor of three relative to the resolved one, and the use of calorimeter 
islands (see next section) instead of cells. For the HERWIG systematic study, a fraction of 
the events was generated with pTmin set to 5.0 GeV. The systematic errors were evaluated by 
calculating the root mean squared deviation between the original data point and each of the 
points found using the systematic variations mentioned above, weighting each deviation by its 
statistical errors. In the following figures the statistical errors are shown as the inner error bars; 
the outer bars show the statistical and systematic uncertainties added in quadrature. 



8 Photon remnant properties 

In Fig. 3, we show the pseudorapidity, transverse momentum (with respect to the Z axis), 
and energy distributions for the photon remnant, corrected back to the hadron level. The 
average correction factors are typically 1.2 and are approximately constant for each of the 
three variables. 

The 773 distribution is shown in Fig. 3a. The corrected data and the expectations from PYTHIA 
(solid histogram) disagree in the negative 7/3 region, as observed in Fig. 2d. The measured 
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distribution peaks at higher values of 773 than the Monte Carlo prediction. A similar effect 
can be observed in the transverse momentum distribution (Fig. 3b), where the requirement 
773 < — f is applied for this and all following figures. Here also, the data show a higher average 
value. The distribution peaks at 1.5 GeV with a tail extending to 6 GeV. The mean value of 
the photon remnant pt at the hadron level is measured to be (PT3) = 2.1 ±0.2 GeV. The Monte 
Carlo expectation is (p T3 ) = 1.44 ± 0.02 GeV. 

The energy distribution (Fig. 3c) peaks around 7 GeV and extends to 20 GeV. The solid 
histogram shows the Monte Carlo expectation which agrees with the data, except in the lowest 
bin where the Monte Carlo expectation is about 30% higher. 

The unlikely possibility has been studied that an incorrect Monte Carlo description of the energy 
response in the RCAL might result in the observed discrepancy between the data and the Monte 
Carlo expectation. The RCAL energy scale was reduced by 10% within 10° of the RCAL beam 
pipe to simulate additional inactive material. The discrepancy in the pT3 distribution was not 
improved, and the disagreement with the data for the £"3 distribution became worse. 
The observation that the average transverse momenta of the photon remnant is higher than 
expected from PYTHIA, is in qualitative agreement with the theoretical predictions mentioned 
in the introduction. Therefore, following [|1(J, we have compared the data with Monte Carlo 
events generated with a harder intrinsic transverse momentum spectrum for the partons in the 
photon, i.e. dN/dtf oc l/(k 2 t + k%). The "PYTHIA high k t " results are shown as the dotted 
histograms in Fig. 3. The parameter ko is determined by minimizing the x 2 between the Monte 
Carlo hadron level and the corrected data pT3 distributions. The result is ko = 0.66 ± 0.22 
GeV. This corresponds to (k t ) ~ 1.7 GeV, as compared to 0.4 GeV for PYTHIA with default 
parameters. The agreement with the energy distribution is unchanged, while it is considerably 
improved for the r] 3 and ptz distributions. On the other hand, the distributions for the first 
two clusters remain unchanged (not shown). Equally good agreement between the data and 
the Monte Carlo simulation has been achieved by reweighting the Monte Carlo events to the 
default intrinsic k t spectrum (not shown). In this case, k = 1.90±0.21 GeV, much higher than 
the default value of k = 0.44 GeV assumed in PYTHIA. These results show that adjusting the 
intrinsic transverse momenta of the partons in the photon is a way to improve the agreement 
between the data and the Monte Carlo predictions. 

One of the well-known properties of jets is that the average energy transverse to the jet axis is 
limited as the jet energy increases. In general, this results in an average energy transverse to the 
jet axis per particle, (Ej), of the order of a few hundred MeV. In measuring this quantity for 
the photon remnant, we used calorimeter islands. Islands are groups of calorimeter cells more 
closely related to particles than individual calorimeter cells and therefore represent a better 
choice for this measurement. The analysis performed with calorimeter cells provides an upper 
limit on the systematic uncertainties of this measurement. Figure 4a shows the average value of 
{Ej) versus the cluster energy for the third cluster, both for the data after corrections and for 
the Monte Carlo simulation. The corrections include the effect of particles lost down the beam 
pipe. The largest component of the systematic errors (+30%) comes from using calorimeter 
cells instead of islands. 

The mean value of (Ej) starts below 200 MeV, and slowly increases with the remnant energy, 
E 3 . The total cluster energy, on the other hand, spans a range from 3 to 21 GeV. This 
result demonstrates that the photon remnant exhibits limited transverse energy per particle, 
as expected for jet-like objects. Further support to this conclusion is given in Fig. 4b, which 
shows the average values of the corrected total transverse (Eji?^) and total longitudinal (EjE^) 
energy of the third cluster, with respect to the cluster function of the energy of the 
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cluster. The longitudinal component increases much faster than the transverse energy, and 
most of the cluster energy is along the cluster axis. This is consistent with a jet-like structure 
of the photon remnant. 

We also studied how the energy is distributed around the axis of the third cluster. Figure 
4c shows the corrected energy flow of the third cluster as a function of 1 — cosO for both 
the data and the Monte Carlo simulation. Here, G is the angle of the particle with respect 
to the cluster axis. This is effectively a plot of the energy deposited in rings of fixed area 
centred on the cluster axis. Because there is not a simple correspondence between particles and 
calorimeter cells, it is difficult to construct a correlation matrix between the generated (hadron 
energy) and experimental (calorimeter cell energy) quantities. Therefore, these distributions 
are corrected bin-by-bin. The average correction factor is around 1.15. In Fig. 4c, we have 
required reconstructed (uncorrected) jet energies between 8 and 14 GeV and hadron level jet 
energies between 8 and 15 GeV. The statistical errors are the error on the mean. The energy 
distribution for the data is quite collimated. The Monte Carlo simulation agrees very well with 
the data, indicating that the fragmentation of the remnant is understood. 



9 Comparison of the photon remnant with jets from 
hard scattering 

Having established the jet-like properties of the photon remnant, we next compare it with the 
jets originating from parton hard scattering. The comparison of these two types of jets is of 
interest because one is the debris of the photon and is a low-pr jet, with p T typically well below 
6 GeV, while the other two jets come from the hard scattering of the partons in the photon 
and proton and are high-pr jets, with a minimum of 6 GeV. A comparison between low-p^ 



and high-px jets has been proposed as a method of studying the hadronization process [29 . 
Figures 4d, e, and f present the results of this comparison for the average value of (Ej,), {Y^E^) 
and (EiE l L ), and 1 — cosO, respectively. In Fig. 4f, we again require reconstructed (uncorrected) 
jet energies between 8 and 14 GeV, and hadron level jet energies between 8 and 15 GeV. This cut 
is especially important for this figure in order to compare jets with approximately equal energies. 
In all figures, good agreement between the third cluster and the two hard jets is observed. From 
these comparisons we conclude that, in the kinematic region and for the variables studied, the 
low-pr photon remnant jet exhibits the same hadronization characteristics as the high-p^ jets 
originating from the hard scattering process. 



10 Conclusions 

For the first time, in a sample of quasi-real photon-proton collisions, the photon remnant 
produced in resolved photon interactions has been isolated. The selected events contain two 
high-p T jets with p T > 6 GeV and r\ < 1.6, and 130 < W 1P < 270 GeV. The properties of the 
photon remnant, as defined by a cluster with i] 3 < — 1 and E 3 > 2 GeV, are studied and shown 
to exhibit a collimated energy flow with a limited transverse energy with respect to the cluster 
axis, characteristic of a jet structure. 

The leading order QCD Monte Carlo simulation, with default parameters, does not reproduce 
the pseudorapidity distribution or the transverse momentum distribution (with respect to the 
incident photon) of the photon remnant. The mean value of Pt for the photon remnant, 
2.1 ±0.2 GeV, is substantially larger than the Monte Carlo expectation. Better agreement can 
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be obtained by increasing the average intrinsic transverse momenta of the partons in the photon 
to about 1.7 GeV. These results are in qualitative agreement with theoretical expectations of 
substantial mean transverse momenta for the photon remnant. 

The photon remnant has also been compared, in the laboratory frame, with the two high-p^ 
jets originating from the parton hard scattering. Although the origins of these two types of jets 
may be quite different, within the present statistics and in the kinematic range studied, they 
exhibit similar properties for the energy flow and the transverse and longitudinal energy with 
respect to the jet axis. 
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Figure 1: Examples of leading order diagrams for (a) direct and (b) resolved hard photopro- 
duction interactions. 
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Figure 2: Pseudorapidity distributions of clusters for the inclusive photoproduction sample. 
The three clusters are sorted by p T cal with (a) having the highest, (b) the second highest, 
and (c) the lowest pT cal ■ Each Monte Carlo distribution is independently normalized to the 
data in the region if al < 1.6. In (c) and (d) the direct contribution alone is shown as the 
dashed line. The rf al distribution of the lowest p T cal cluster (as in (c)) is shown again in (d) 
after requiring prTi > 5 GeV, rf^ < 1.6, and E^ al > 2 GeV. Resolved events are selected by 
requiring r]^ al < — 1 (vertical line in (d)). The dotted line shows the expectation using the 
LAC1 photon parton parameterization. 
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Figure 3: (a) Pseudorapidity distribution of the third cluster corrected back to the hadron 
level, (b) Corrected distribution, (c) Corrected energy distribution. In (b) and (c) we require 
i] 3 < — 1. The solid histograms are the hadron level distributions given by the default version 
of PYTHIA (Gaussian with k = 0.44 GeV). For each of the figures, the dotted line shows 
the Monte Carlo predictions with dN/dkf oc l/(/cf + k$) and k = 0.66 GeV, corresponding to 
(kt) ~ 1.7 GeV. 
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Figure 4: (a-c) Comparison between the data and the hadron level Monte Carlo expectations, 
(a) The mean value of (Ej,), the average energy transverse to the cluster axis per particle, as a 
function of the cluster energy, (b) The average values of the total transverse (Eji?^) and total 
longitudinal (EjE^) energy, (c) The flow of energy around the cluster axis, (d-f) Comparison 
between the photon remnant (cluster 3) and the two high-p^ jets. The error bars show the 
systematic and statistical errors added in quadrature. 
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